Polaronic memristor strongly coupled to electrodes 
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j*— «j , Attractive electron correlations due to an electron-vibron interaction (EVI) can overcome the 

f—*) . direct Coulomb repulsion of polarons in strongly deformable molecular quantum dots (MQDs). If 

' it realizes, a switching appears in the I-V characteristics of the degenerate nonadiabatic molecular 

bridges weakly coupled to electrodes providing a route to ultrafast 'memristors' (memory-resistors) 
as the basis for future oscillators, amplifiers, and other important circuit elements. Here, we extend 
our theory of polaronic memristors to adiabatic MQDs strongly coupled to the leads and show 
that the degeneracy of MQD (or other multilevel energy structure) along with the polaron-polaron 
attraction is a necessary ingredient of its switching behavior in the strong-coupling regime as well. 
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PACS numbers: 63.22.Gh, 62.23.Hj, 73.23.Hk, 71.38.Mx 



Different nano-size devices are being proposed and investigated [H, 0J3| that exhibit some kind of current 'switch- 
ing' behavior 

am, 

'negative differential resistance', and 'memory' [fj. There is currently a surge of interest in 
various systems showing memristor behavior (see, for example [3] and references therein) that can potentially be used 
for Resistive Random Access Memories (RRAMs). The practical significance of determining a precise microscopic 
mechanism of such a behavior is difficult to overestimate. 

Polarons - electrons strongly coupled with lattice vibrations (phonons or vibrons) - play a key role in the transport 
and optical properties of many systems of reduced dimensions and dimensionality [2 [t| [l(| • They may provide an 
almost instantaneous switching mechanism [ll[ , when a bistable current state of vibrating nano-circuits appears due 
to the attractive electron-electron correlations ('anti- Coulomb blockade' of the insulating state [HI, [H|), if MQD 
is many-fold degenerate. A mean-field approximation with respect to EVI [lij fails to describe intrinsic molecular 
switching, so that the exact analytical or numerical solutions of a polaronic transport problem are required [III Il5j. 

An exact analytical solution of the problem fully accounting for the strong EVI and polaron-polaron correlations 
in MQD has been found by us in a weak molecular-lead coupling limit, where an inverse lifetime, Y/h, of a polaron 
on the dot is small compared with the characteristic vibron frequency, Y/h <C coq [ll| (nonadiabatic regime). Here 
we extend our theory of polaronic memristors to MQDs strongly coupled with the leads, where Y/h ^> wo (adiabatic 
regime) . 

The molecular Hamiltonian includes the Coulomb repulsion, U c ' , and the electron-vibron interaction as [H| | 



> 

o 

m 

(N ■ ^ i ^ 

+ Yl w «M7«d9 + H.c.) +J2^i(4d q + 1/2). (1) 



Here e M are energy levels of a rigid molecule, = c^c^ is the molecular occupation number operator, d q annihilates 
vibron, to q is the vibron frequency, and ^ m is the EVl constant (q enumerates the vibron modes, and we have taken 
h = ks = 1 here and below). This Hamiltonian conserves the occupation numbers of molecular states n M . When the 
molecule is attached to two (for simplicity symmetric) leads, the complete Hamiltonian H = H m +Hi+J2^ k *fc/i c /*( a fc+ 

bk) + H.c. includes the hopping to the leads, and the Hamiltonian of the leads Hi = J^fc£fe( a fc a fc + &!&fc)> where au 
and bk annihilate the electron in the left and right electrodes, respectively. 

With some conventional assumptions (energy independent width Y = 27r^ fc f^ <5(i£ — and quasi-equilibrium 
electron distributions in the leads ) the current through MQD is given by the Landauer-type expression (or by the 
Fermi golden rule) as [111 ]: 

POO 

I(V) - I / dE [h{E) - f 2 (E)} p(E), (2) 



allowing for a transparent analysis of the essential physics of the switching phenomenon. Here Iq = eY, fi(2)(E) = 
l/[exp[/3(£ , + A=FeV r /2)/T] + l] is the electron distribution function in left (1) and right (2) metallic leads, respectively, 
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FIG. 1: Current- voltage characteristic of the nondegenerate (d — 1) MQD weakly coupled to the leads at T = OK, Wo/A = 0.2, 
and j 2 — 11/13. The vibron ladder in clearly seen in the I-V characteristic, but there is no hysteresis. 



[3 = l/T the inverse temperature, A the position of the lowest unoccupied molecular level with respect to the Fermi 
level at V = 0, and V the voltage drop across MQD. 

The molecular density-of-states (DOS) p(E) = — (1/tt) ^2 ImGj} '(E) generally depends on all interactions and 

molecule-lead hopping integrals as determined by the Fourier component G^(E) of the retarded Green's function 
(GF) G^(t) = —iQ(t) (c^(t)c|, + c^c^t)) , where c M (i) is the Heisenberg annihilation operator. For calculating p(E), 
one can neglect molecule-lead coupling in the weak-coupling regime while keeping all orders of the Coulomb repulsion 
and EVI by means of the canonical displacement transformation of the molecular Hamiltonian |T]). In particular, 
assuming EVI with a single vibronic mode, oj q — loq and the coupling ^ m = 7, and constant Coulomb integrals, 
[/?,, = V c , one obtains for a <i-fold degenerate single-level dot (Tlj . 
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e ^r- [(1 _ n )S(E - rU - Ilu ) + nd(E - rU + lu )] 
+ (1 - 5 m )e-^[n5{E - rU - Iuj ) + (1 - n)S(E - rU + Ilo q )} 
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Here Z = exp [— I7I 2 coth(/3fio;o/2)] accounts for a familiar polaronic renormalization of the hopping integrals, £ = 
7I 2 / sinh(/3fiwo/2), // (£) is the modified Bessel function, and 5ik is the Kroneker symbol (/, k = 0, 1,2, ...), and we 
have taken the position of the level as zero, = 0. The resulting polaron-polaron interaction, U = V c — 2E p , where 
E p = Ij^ujq is the polaron level shift, comprises the Coulomb repulsion, V c , and the vibron-mediated attraction. 
An important feature of DOS ([3]) is its nonlinear dependence on the occupation number n of degenerate molecular 
states owing to the correlation side-bands with the spectral weight Z r (n) = (ef— l)!n r (l — n) d ~ 1 ~ r / (r\(d — 1 — r)!). 
The DOS, Eq.([3]), contains full information about all possible correlations in transport, in particular, the vibron and 
correlation side-bands. It is derived by solving the finite system of coupled equations for A^-particle Green's functions, 
as described in Ref. [TT[ | . 

Equating incoming and outgoing numbers of electrons in MQD per unit time, one obtains the self-consistent equation 
for the molecular state occupation number n as [ll| : 

2nd = J dEp (E) [fx (E) + f 2 (E)} , (4) 

which automatically satisfies the condition ^ n $J 1. In particular, for the nondegenerate MQD (d = 1) and T = OK 
the result is 
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where ojf = ZY,Zoi\l\ l l W{^0 - A + eV/2) ± 0(lwo - A - eV/2)], and 6jf = £E»=o(l7r/J0[6(-k*> - A + 
eV/2)±©(— Iujo — A — eV/2)]. Here 0(:c) = 1 if x > and zero otherwise. The current is single- valued, Fig.[TJ with the 
familiar ladder due to the phonon side-bands. The double-degenerate level d = 2 obviously provides greater number 
of elementary processes for conductance reflected in a larger number of steps on the phonon ladder in comparison with 
d = 1 case. However, the current remains single- valued showing no switching behavior [ToL ITlj in the weak molecular- 
lead coupling regime, contrary to the mean- field approximation (MFA) of Ref. [3], which yields an artificial switching 
behavior for a non-degenerate level, d = 1, in this regime. In fact, the switching appears only in the I-V curves of 
many- fold degenerate dot when U < and the temperature is low enough [l2|. 

We analyze below an effect of a strong molecular-lead coupling r > wo on these results. When EVI is also large 
(7 3> 1 ) , it could not be treated perturbatively, so that one has to sum all diagrams of the perturbation expansion in 
powers of 7 in order to calculate GFs. As suggested in Ref. (lr| . this can be done with a lowest-order non-equilibrium 
linked cluster expansion (NLCE) [13, which provides an approximate resummation of the whole series for GFs as: 



G(t) = J2^W r (t) « G Wexp[ 7 2 F(i)]. 



(7) 



Here Go(t) should be calculated without EVI [i.e. with j M = in Eq.([T])], while F(t) is found using the lowest 
(second)-order diagrams F(t) = Gq 1 {t)W\{t) . 

Calculating zero Go (t) and second order (in 7) W\ (t) is a textbook exercise for the non-degenerate molecular level. 
Using equations of motion, ic — t^2 k (a k + fyfc) a- n d icik = + tc, ib k = Ck^k + tc, one can readily find the Fourier 
component of the retarded GF of the rigid non-degenerate dot as G§(E) = (E — iT)^ 1 , and a lorentzian zero-order 
molecular DOS, 
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Substituting Eq.® into Eq.f2]) yields the current through a rigid bridge as 
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shown in Fig. [2ji at T = OK for a few inverse lifetimes. 

Using the zero-order lesser, G^(t) ms in exp(— T\t\), and greater, Gq({) w — i(l — no) exp(— T\t\), GFs Jl6| one can 
readily calculate the second order contribution Wi(t). Here uq = [/i(0) + /2(C))] /2 is the zero-order population of 
the molecular level. Then applying NLCE, Eq.([7]), one gets the retarded GF for the strong molecular-lead coupling, 
r ^> ojo, summing infinite number of EVI diagrams (see for details Ref. [H, G R (E) — (E + 2noE p — ir)^ 1 , and 
the molecular DOS, 



P(E) 



Using Eq. (fT0|) in Eqs.fjl]),©, one easily obtains: 
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for the molecular-state population and 
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(10) 
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for the current, respectively, with A = A — 2uqE p = A — E p Q(eV/2 — A) at low temperatures, T <C T. The level 
population, Eq. pTj) . of the non-degenerate dot strongly coupled to the leads, and the current, Eq. (fT2")) , Fig. [2b, remain 
single- valued showing neither switching nor negative differential resistance similar to the weak molecular-lead coupling 
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FIG. 2: Current- voltage characteristics of the nondegenerate (d = 1) rigid molecular bridge at T = OK (a), and of the non- 
degenerate adiabatic (u>o <C T) molecular bridge with the polaron level shift E v — A/2 (6) for different couplings to the leads, 
T/A = 0.01 (the steepest curve), 0.1, 0.2, and 0.5. 



regime, Fig. [T] Taking into account the coupling with the leads beyond the second order washes out the vibron ladder 
in the I-V characteristic, so that the EVI affects the FV characteristics only marginally when compared with the rigid 
dot, Fig. in the strong-coupling (adiabatic) regime. 



As we have shown earlier 15[ , the mean-field approximation (MFA) of Ref. [14j replacing the occupation number 



operator fi in EVI by the average population n erroneously leads to a non-linear equation for n and a spurious 
switching of the non-degenerate dot, which has no physical meaning. The authors of Ref. [1 61 ] have confirmed our 
conclusion for the weak molecular-lead coupling, but argued by using NFCE that, surprisingly, their MFA works well 
in the strong coupling regime, so that the non-degenerate dot is multi-stable contrary to our present result, Fig.[2jD. 
The discrepancy originates in an erroneous replacement of the zero-order occupation no in Gq ,:> (i), G R (E), and as the 
result in Eq. (|10[) . by the exact n leading to double-counting of EVI diagrams in Ref. [16]. In fact, multiple solutions 
for the steady-state population of the non-degenerate dot found in Ref. [13] and Ref.[l6j are artefacts of MFA and 
mis-application of NLCE, respectively. The physical mechanism of the switching is provided by an attraction of two 
polarons on a multiply-degenerate MQD [0) [HI G3] > which is missing in the non-degenerate 'polaron model' of Refs. 
[Til [r| due to the Pauli exclusion principle. 

Treating the polaron-polaron correlations in real multi- level bridges and the molecular-lead coupling at the same level 
of approximation is a challenging problem, which might require numerical techniques. For instance, the conductance of 
deformable molecules has been numerically studied in the framework of a two-impurity Anderson model with positive 
and negative electron-electron interactions and in the two-impurity Anderson-Holstein model with a single phonon 
mode [181 ] . where the spin and charge Kondo effects can occur simultaneously at any coupling strength. Since the 
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FIG. 3: Steady-state population n and the current I /Io of the many-fold degenerate negative— U dot at T — OK with different 
values of the coupling to the leads, and the effective attraction U = — A/2. 



vibron ladders in I-V curves are washed out by the strong coupling with the leads, Fig. [2p, we consider here an effect 
of the molecular-lead coupling on the switching in the many-fold degenerate (d 3> 1) negative Hubbard U model of 
Ref. [T2| allowing for a simple analytical solution, 



Hm = \ U Yl K^' +J2tk(a{a k + b{b k ) +J2 t kv( a k 



H.c. 



(13) 



For the many- fold degenerate dot one can approximate the exact two-body interaction in the Hamiltonian Eq.([T 



by a single-particle self-consistent Hartree-type potential as ^^n^n^i m UY] 



Then, using equations of motion yields the DOS as in the non-degenerate rigid dot, Eq.((8]), but shifted in energy by 
the attraction potential, 
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From Eq.(j4]), the molecular state population for T <C V is: 
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where U = U(d- l)/2. 
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For r = 0, there are two stable solutions of Eq. (fT5|) , n = and n — 0.5 in the voltage region of bi-stability, 
A - \U\ < eV/2 < A, and only one solution, n = 0.5 for eV/2 > A, and one solution n — for eV/2 < A - \U\. The 
coupling with the leads does not destroy the many-fold degenerate memristor as long as T is less than some critical 
value, r < r c , which is small compared with A, Fig. [3J When T -C A, one can replace the last arctan in Eq. (fT5]) by 
7r/2 and resolve this transcendental equation with respect to the bias voltage as eV/2 « A + 2Un — V cot(27rn). The 
boundary of bi-stability is found from dV/dn = , which yields the terminal value of the level width T c — \U\/tt. The 
dot population and the current, found from Eq.@ using Eqs. (|14ll5[) remain double-valued in the voltage region of 
bi-stability, Fig. [3J where both high-current and low-current branches are stable, while the intermediate branch with 
dn/dV < 0, Fig.[3ji, is unstable [12]. Importantly, the coupling narrows the voltage range of the hysteresis loop, but 
the transition from the low (high)-current branch to the high (low)-current branch remains discontinuous as long as 
r < r c . The Hartree-type approximation of the correlation potential is fully justified by the exact solution of Ref. 
[l2T ] in the weak molecular- lead coupling regime. 

The actual mechanisms of nanocircuit switching are of the highest experimental and theoretical value. Further 
progress will depend upon understanding of intrinsic mechanisms of nano-memristors. Here, we have extended our 
multi-polaron theory of the current controlled switching mechanism to the adiabatic molecular-size bridges with a 
significant coupling to the leads (r>wo). The degenerate (or multilevel) MQD shows the hysteretic memory, if the 
degeneracy of the level (or the number of levels) is large enough, the Coulomb repulsion is over-screened by EVI or 
by any other attractive mechanism, and the coupling to the leads is below the critical value, Fig. [31 Importantly, the 
switching does not exist in non-degenerate molecular dots neither weakly, Fig. 1, nor strongly coupled to the leads, 
Fig- (2b, contrary to the claims in Refs. [1, [lj|, [l6|. In all cases, bi-stability is destroyed by a very strong coupling to 
the leads, T > r c . 
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